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ABSTRACT 

We perform population synthesis modelling of a magnitude-limited sample of 4391 
Luminous Red Galaxies selected from the Sloan Digital Sky Survey Data Release 4 
(SDSS DR4). We fit measured spectral indices using a large library of high resolution 
^ ' spectra, covering a wide range of metallicities and assuming an exponentially decaying 

I star-formation rate punctuated by bursts, to obtain median-likelihood estimates for 

lO ■ the light-weighted age, metallicity, stellar mass and extinction for the galaxies. The 

' ages lie predominantly in the range 4-10 Gyr, peaking near 6 Gyr, with metallicities 

in the range -0.4 < [Z/H] < 0.4, peaking at [Z/H] w 0.2. Only a few per cent of the 
spectra are better fit allowing for a burst in addition to continuous star-formation. 
For these systems, typically one-quarter to one-third of the stars are formed in the 
burst. The total stellar masses of all the galaxies are confined to a very narrow range 
around ~ 3 x 10^^ Mq, consistent with the expected homogeneity of the sample. Our 
Qh' results broadly agree with those of previous groups using an independent population 

Q . synthesis code. We find, however, that our choice in priors results in ages 1-2 Gyr 

' smaller, decreasing the peak formation epoch from about z = 2.3 to z = 1.3 for 

^ I the stars. To describe the distribution in measured mean metallicity of the galaxies, 

. we develop a metal evolution model incorporating stochastic star-formation quenching 

motivated by recent attempts to account for the apparent 'anti-hierarchical' formation 
of elliptical galaxies. Two scenarios emerge, a closed box with an effective stellar yield 
' of 0.26, and an accreting box with an effective stellar yield of 0.10. Both scenarios 

5-H ' require an IMF weighted towards massive stars. They also require characteristic star- 

formation quenching times of about 10® yr, the expected lifetime of luminous Quasi- 
Stellar Objects. The models predict an anti-correlation between the age and mean 
metallicity of the galaxies similar to that observed. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters - galaxies: stellar content 



1 INTRODUCTION 

In the hierarchical merger scenario of structure formation, 
more massive systems form from the mer ger of less massive 
ones (|Peebles fc Dickelll968l : |Peebleslll983 ). While numerical 
simulations of the evolution of dark matter halo es in a Cold 
Dark Matter dominated universe bear this out (jPavis et al.l 
1 19851 ). predictions in the hierarchical scenario for galaxies 
are co mplicated by the nee d to account for gaseous dissi- 
pation l| White fc Reesll 19781 ). without which galaxies would 
lose their integrity at each level of the merger hierarchy. 
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The most straightforward prediction of the hierarchical sce- 
nario is that smaller mass galaxies form prior to the more 
massive, although the effects of galaxy mergers, as distinct 
from halo mergers, cou ld substantially complicate this pic- 
ture dBlumenthal et al][l984i ). Early semi-analytic estimates 
of the formation and properties of galaxies have suggested 
that half of the stars in the Universe formed at r edshifts 
z < 1.5, and only som e 10 per cent at 2; > 3 (Co le et"ai] 
ll994l : lBaugh et aflligg^ ) . Elliptical galaxies in particular are 
predicted to show a trend of decreasing age with increas- 
ing luminosity, w i th the brightest ellipticals the youngest 
jKauffmannll 19961 : iKauffmann fc Chariot! Il998l ) . 

Observations have generally not supported the pre- 
dicted ages and their trends with galaxy properties. The 
contrary evidence is derived primarily from two methods: 
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measurements of the luminosity function of galaxies and 
its evo lution, and age esti mates based on population syn- 
thesis. ICowie et al ] nig?) discovered that the character- 
istic luminosity of galaxies undergoing rapid star forma- 
tion has been declining with t ime at redshif t s z < 1.7, a 
result extended to 2 « 3 by ISawicki et aP (|l997h . As a 
function of galaxy class, the later the t ype of galaxy, th e 
later is its peak period of star formation (|Hevl et al.lll997t ). 
Since then, an abundance of evidence for anti-hierarchical 
galaxy formation (or so-called 'downsizing') has accrued 
based on galaxy number counts and colours, generally indi- 
cating that the 'quenching mass', at which the major episode 
of st ar formation ceases , decreases with decreasing red- 
shift (iPozzetti et al.ll2003l: lFontana et al.ll2004|: iBundv et al 



2005| : iDrorv et all l2005riBundv et all I2OO5I : Tcimatti et al 
20061 ). Massive galaxies appear to undergo a higher merger 
rate at high redshifts ( z » 3) than a t low redshifts, set- 
tling down by 2 « 1.5 (|Conselicell200i ). again in apparent 
opposition to the hierarchical merger scenario. 

The picture that appears to be developing is that star 
formation occurred early in systems with a wide range of 
masses, but was terminated through some unknown process, 
first in the progenitors of the most massive galaxies, and 
progressing steadily towards less massive systems with time. 
The progenitors of massive galaxies (with stellar masses to- 
day in excess of 10^" M© ), subsequently merged, possibly 
with a small amount of additional star formation, but set- 
tled into largely quiescent systems by z — 1.5. 

Numerical modellers have sought to catch up with these 
findings, drawing increasingly on the role of feedback. While 
supernova feedback has long been assigned a key role in 
regulating star formation, Active Galactic Nuclei have now 
been invoked as a means of terminating star formation in 
the most massive systems, and proposed as the origin of the 
dichotomy between blue galaxies and the more massive red 
iDekel fc Birnboim|[2006h . Incorporating feedback by super- 
novae and Active Galactic Nuclei has become a promissing 
means of reconciling the mass-sensitive evolution in num- 
ber coun ts, the m erger rate s, and the ages of stel lar pop- 
ulations (|Bower et al.. .20061 ; iLucia fc Blaizo^ bOOej ). These 
models may be expected to develop further as knowledge 
of the star formation history of galaxies of different masses 
and in different environments continues to grow increasingly 
precise. 

In the past, population synthesis model estimates of 
ages were long plagued by the age-metallicity degeneracy. 
Improvements in the stellar atmosphere models, the data, 
and the understanding of the sensitivity of absorption lines 
to age and metallicity have since led to more relia ble age and 
metallicity determinations. iTrager et all (|2000bm found a 
range of ages for elliptical galaxies in the field, with lower 
velocity-dispersion galaxies tending to be younger. By con- 
trast, ellipticals in clusters tend to be old. They also found 
that larger velocity-dispersion galaxies ten d to show in- 
crease d suppression of iron-peak elements. ICaldwell et al.l 
(|2003l ) similarly found that smaller velocity dispersion sys- 
tems tend to have younger luminosity-weighted ages. Evi- 
de nce for rece n t bur st activity in field ellip ticals was found 
bv lTreu et all (|2002l ). iThomas et al.l (|2005l ) have suggested 
a shift in the peak star formation epochs of early type galax- 
ies, with massive early-type galaxies in low-density environ- 
ments ~ 2 Gyr younger than their high-density environment 



counterparts, and the star formation rate peak shifting to 
lower redshifts the lower the mass of the ga laxy. Old ages 
for cluster ellipticals were similarly found bv lYamada et all 
(2006). 

Estimates of the current to mean star formation rates, 
or 'birthrate,' of galaxies using population synthesis mod- 
els for galaxies of a range in mass and type show a 
tren d of decreasing birthrate with increasing stellar mass 
(.BoseUi et al.l l200l[ iGayazzi et al.] |2002|: iKauffmann et all 



I2OO3I : iBrinchmann et al.ll2004l : iGallazzi et all 120051 ) . again 
indicating that the more massive a galaxy is, the earlier 
most of its stars were formed. 

The confluence of wide-field surveys and high-speed 
computing is now making feasible the testing of galaxy 
formation scenarios by detailed population synthesis mod- 
els with a spectral resolution and coverage previously im- 
possible on a large scale. In this paper, we apply popu- 
lation synthesis models to estimate the mass, metallicites, 
and star formation histories of Luminous Red Galaxies 
LRGs) drawn from the Sloan Digital Sky Survey (SDSS) 
York fc et alll2000l ). We have several reasons for concen- 
trating on LRGs: 1. they form a large sample of galaxies 
chosen according to well-defined selection crit eria forming a 
well-defined sample (jEisenstein fc et al.ll200ll ). 2. the galax- 
ies form a uniform population both photometrically (|Fabeii 
|l973; Visvanat han fc Sandagelll977l). and spectroscopically 
(fpaber .1973 : .Efisenstein fc et ahF 20031 ). 3. the galaxies are 
bright, providing high signal-to-noise ratio spectra, 4. they 
represent the extreme high mass end of galaxy stellar masses 
and an extreme in anti-hierarchical behaviour, 5. star for- 
mation has ceased in all but a few percent of the systems 
(jRoseboom fc et al.ll200(} ). and so the influence of emission 
lines associated with ongoing star formation on the key ab- 
sorption features used to infer the properties of the galaxies 
is minimal, 6. they are an intrinsically interesting class of ob- 
jects, being the most massive and most strongly clustered 
galaxies in the Universe. 

Population synthesis analyses for galaxies in SDSS 
Data Release 2 (DR2) and ear li er da ta were performed 
previouslv bv [Kauffmann et all ||2003|) (he r eafter K03), 
IBrinchmann et al.l (|2004l ) and IGallazzi et al.1 l|2005h (here- 
after G05). These studies provide useful comparisons for 
our results. Our analysis, however, differs in several re- 
spects from these earlier efforts. We exploit high resolution 
synthetic spectra to model several absorption indices that 
match the resolution of th e spectra. The analyses of K03 and 
IBrinchmann et al.l l|2004) used the 4000- A break (D„4000) 
and Balmer absorption, or emission lines, to estimate the 
star formation rates and histories. The closest population 
synthesis analysis to ours is that of G05, who similarly based 
their model-fitting on several high resolution absorption fea- 
tures. This work, however, differs from theirs as follows: 1. an 
alternative population synthesis model is used, utilising an 
independent set of stellar templates to construct the spectra 
and a different set of stellar evolution tracks, 2. the sample 
of LRGs analysed is substantially larger, and 3. we use a 
different set of absorption features specifically adapted to 
the LRG sample analysed rather than the full galaxy pop- 
ulation analysed by G05. In consequence, our results are 
complementary to theirs, and provide independent checks 
on their findings, while at the same time allowing us to ex- 
tend the results using the larger sample. G05 also found 
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some discrepant absorption indices poorly matched by the 
best-fitting models. Our alternative set of models sheds light 
on the origin of these discrepancies. 

Another goal is to quantify the metallicities of the 
systems and their deviations from solar abundances. 
Previous studies have suggested that gia nt ellipticals 
and LRGs contain super-solar m etallicities ijFabeij [19731 : 
IPickled [l985l : iGallazzi et al]l2005l). with indications of a- 
enhanced abundances jWorthev etall 1 19921: iDavies et al.l 
1 19931 : iTraeer et all l2000bl lat lEisenstein fc et all l2003h.~a r 
though agreement is not universal I Kelson et al.l 2006t ). 
The galaxies also show a trend of increasing metallic- 



ity with veloci t y dispersion or mass (ITraeer e t al. 2000a" 
Caldwell et aD l2003l : I Thomas et al.l 



20051 : iKelson et al 



20061 ). Little attention has been devoted to the modelling 
of the metallicity distribution of early type galaxies. In 
Section |4l we present models incorporating into standard 
closed box and accreting box models the sudden quenching 
of star-formation, as may occur if Active Galactic Nuclei 
heat the gas, effectively shutting down gas cooling, to de- 
scribe the evolution of metallicity. We find that the mea- 
sured metallicity distribution provides added support to the 
star-formation quenching scenario in the context of these 
enrichment models. 

In the next section, we discuss the sample selected. In 
Section 3, we describe our methods. Results from the anal- 
ysis are presented in Section 4. A discussion of the results 
and conclusions are provided in Section 5. 



2 SAMPLE SELECTION 

The spectra used in this analysis are taken from 
the Sloan Digital Sky Survey D ata Release 4 (DR4) 
(jAdelman-McCarthv fc et al.l l2006h . The SDSS is a five- 
band imaging survey of the northern Galactic gap, com- 
bined with an extensive spectroscopic follow-up program. 
The spectra cover a large wavelength range (3800-9200 A) at 
high resolution A/ A A ~ 1800 taken with 3 arcsec diameter 
fibres. The target spectra are selected using the Luminous 
Red G alaxy (LRG) cuts as defined in lEisenstein fc et aP 
()200ll ). The aim of the SDSS LRG survey was to produce a 
'volume-limited' sample of intrinsically luminous (L J5 31/*), 
intrinsically red galaxies out to redshifts of z — 0.5. Such a 
sample is useful as a probe of large scale structure in the 
Universe. Here we use the data to study the properties of 
giant elliptical galaxies. 

The LRGs are selected from the SDSS database by re- 
quiring the TARGET_GALAXY_RED flag to be set. An additional 
cut requires the redshift to fall within the range 0.15 < z < 
0.4, with the lower limit in place to ensure that underlumi- 
nous galaxies do not enter the sample. A velocity dispersion 
cut is also made, requiring 70 kms~^ < a < 450 kms~^. 
This initial selection yields a dataset of 24,615 galaxies. 

Before detailed analysis was undertaken, each spectrum 
was blueshifted to its respective rest-frame and resampled 
in steps of 2 A to match the resolution of the model spectra. 
Unlike in DRl, the DR4 spectra are not pre-corrected for 
foreground Galactic d ust extinction. We ap plied this correc- 
tion using the maps o flSchleeel et al.l (Il998h . combined with 
an extinction law of ICardelli et al.ri|l989l ) and Rv = 3.1. 



All magnitudes used are also first corrected for foreground 
extinction, using values given in the SDSS database. 

Our analysis is based predominantly on fits to five spe- 
cific absorption indices, defined below in Section [3.21 LRG 
spectra from the SDSS, however, contain not only absorp- 
tion features but also emission lines produced by hot, in- 
terstellar gas. These lines can contaminate the stellar ab- 
sorption lines and hence suppress the measured absorption 
index values. Unlike G05, who estimate the effect of emis- 
sion lines by fitting a linear combination of template spec- 
tra, we choose to simply exclude those galaxies where emis- 
sion lines are detected within the index passbands. This 
is done by masking pixels from the index calculation for 
which the SP_MASK_EMLINE flag is set. We also reject pix- 
els which have no data present, no sky information, or 
bright sky levels (flags SP_MASK_NODATA, SP_MASK_NDSKY and 
SP_MASK_BRIGHTSKY). After this, an index wiU only have a 
measured value if there is more than one unmasked pixel 
in each of the three defining passbands. To make the final 
cut, spectra must then have measured values for all five of 
the matching indices described in Section 13.21 These cuts 
reduce our sample to 4391 spectra in total. Selecting this 
subset of galaxies does not systematically affect our results; 
when parameter estimates are made using a subset of the 
five fitting indices, the overall distribution of that param- 
eter is the same as when all five indices are present. This 
suggests our findings for the LRG population as a whole are 
not much affected by contamination from galaxies with ap- 
preciable emission lines. This is expected from the low (few 
per cent) contributio n of actively star-formin g galaxies to 
the LRG population (|Roseboom fc et al.ll2006l ). 



3 METHOD 

3.1 PEGASE synthetic spectra 

To analyse the physical properties of the LRG sam- 
ple, we use the PEGASE evolutionary synthesis code 



(|Fioc fc Rocca-Volmc rangc 19971), version 2.0, coupled with 
a new li brary of high resolut i on Ku rucz stellar spectra gener- 
ated bv I Murphy fc Meik sin' ('2004'). With this version of the 
code, developed bv IMu' dIiv (2003), we can produce spec- 
tra of stellar populations with a resolution of 2 A in the 
wavelength range 3000-10000 A, ideal for comparing with 
the SDSS spectra. The exceptions are for a subset of cool 
stars, for which only 20 A resolution synthetic spectra are 
available, and similarly for very hot stars, although these 
tend to be nearly featureless so that high resolution is not 
cr ucial. Details of the synth etic s pectra used ar e provided 
bv lMurphv fc MeiksinI (|200i) and lMurph^ (|2003| ). 

The PEGASE code uses the standard PAD QVA stel- 
lar evolutionary tra cks of JB rcssan et al.l (119931) . supple- 
mented by those of lAlthaus fc Benvenutd (|l997l ) for low 
mass stars becoming helium white dwarfs. Pseudo-tracks of 
thermally pulsing asy mptotic giant branch (AGB) star s are 
constructed following iGroenewegen fc de Jon j (|l993l ). Fi- 
nally, post- AGB and CO white dw a rf tra cks are taken from 
iBloeckej (| 19951 ) and ISchoenberneJ (|l983l ). The positions in 
the HR diagram of unevolv e d low-mass stars are taken from 
IChabrier fc Baraffd (|l997l ). iBruzual fc Chariot! (2003) find 
that the PEGASE predictions for colours and mass-to-light 
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ratios agree closely with those of their standard population 
synthesis model for simple stellar populations (SSPs) older 
than 1 Gyr when the same IMF is adopted. 

We note that the evolutionary tracks are computed with 
solar abundance ratios, while our galaxies are expected to 
have ratios that differ from the solar values. This will dic- 
tate the choice of stellar diagnostics used to constrain LRG 
properties, and is discussed further in the next section. 

We adopt a Salpeter initial mass function (IMF) for 
masses 0.1 — 120 Mq. T his choice diffe rs fr om those of K03 
and G05, who adopt the lKroupal (|200ll ) and lChabrieil (|2003l ') 
IMFs, respectively. (Several IMFs are provided by Chabrier; 
the form adopted by Gallazzi et al. was based on Hippar- 
cos measurements of disc stars in the Milky Way.) Both 
.roups perform their mod elling using the Bruzual-Charlot 
Bruzual fc CharlotI bo ps') population synthesis code. Tests 
reported in lBruzual fc C hariot (2003) on simple stellar pop- 
ulations show that the Salpeter IMF predicts a mass-to-light 
ratio a factor 1.5-2 greater than these other choices, which 
flatten at the low mass end compared with the Salpeter IMF. 
For masses above a solar mass, both IMFs are very similar 
in shape to the Salpeter IMF. Since the light of the LRGs is 
dominated by more massive stars near the Main Sequence 
turn-off, we expect that the inferred properties of ages and 
metallicities will be largely unaffected by the choice of IMF. 
This is support ed by the similar co lours produced by the 
different IMFs l|Bruzual fc CharlotI |2003). Since there are 
no direct measurements of the IMF in LRGs, we adopt the 
Salpeter IMF for the sake of its widespread usage. 

We also differ from K03 and G05 in the stellar spectral 
libraries used. While we rely on synthesised stellar spectra 
to provide near uniform coverage of the space in surface 
temperature, surface gravity and metallicity, the Bruzual- 
Charlot code incorporates the measured spectra of stars 
from the STELIB stellar spectral library (|Le Borgne et all 
I2OO3! ). The STELIB library has the advantage of having been 
drawn from actual stellar spectra, but at somewhat lower 
resolution (~ 3A FWHM). The coverage in metallicity, sur- 
face temperat ure and surface gravity is also not uniform 
iBruzual fc Ch ariot 2003; Murphy 2003). Our resuHs may 
be viewed as complementary to these earlier efforts, in that 
they test the robustness of the inferred properties of the 
galaxies to the assumptions made in the population synthe- 
sis modelling. We compare our results with those of K03 and 
G05 and comment on any differences in Section [S] Here we 
note that restricting the modelling parameters and galaxy 
sample to those of K03 and G05 results in good agreement 
with their findings. 



3.2 Stellar diagnostics 

In this study we use absorption indices as stellar diagnos- 
tics to model the strength of absorption features in the 
LRG sp ectra. The orig i nal s e t of 21 Lick indices is con- 
sidered (|Burstein et all Il984l : IWorthev et"aLl Il994l '). aug- 
mented with the four high er-order Balmer line indices of 
IWorthev fc Ottaviani' (l99'f) and the high-resolution Hjhr 
index of Jones & Worthcy (1995). Each index is measured 
by calculating the equivalent width over a central bandpass, 
normalised by a pseudo-continuum which itself is defined by 
two adjacent bandpasses. Previous studies have shown that 
some indices are sensitive mainly to age, such as the Balmer 



lines, while others are primarily sensitive to the metallic- 
ity of the stellar pop ulation, such as the Fe and Mg indices 
(|Trager et al.l l2000bl l . E rrors on the i ndices are calculated 
following the method of ICardiel et al.l (1998). 

In addition, we consider the 4000 A break in- 
dex D(4000), which gives the ratio of the average flux 
density in the ba r ids 4 050-4250 A and 3750-3950 A 
(|Bruzual fc CharlotI [ioo^ Also calculated is the narrow- 
band version D„ (4000'l of fealogh et all (1999. ) (4000-4100 A 
and 3850-3950 A), which is believed to be considerably less 
sensitive to reddening effects. The ratio of the strength of 
the Call H and K lines, as integrated fro m 3921-3946 A and 
3956-3981 A, is also considered briefly (jStoughton fc et al.l 
[2002iy Unlike the origin al Lick index studies , which use ana- 
lytical fitting functions (IWorthev et al.iri994l ') to model index 
trends, we follow the same method as G05, directly com- 
paring indices measured from our high resolution synthetic 
spectra to those from the SDSS. 

The stellar spectra used are all based on solar metal 
abundance ratios. External galaxies, as well as having a 
range of metallicities, are expected to have varying abun- 
dance ratios of a elements to Fe-peak elements. It is believed 
that a elements, which include N, O, Mg, Ca, Na, Ne, S, Si 
and Ti, are produced predominantly in Type II supernovae, 
while Fe-peak elements, which include Cr, Mn, Fe, Co, Ni, 
Cu and Zn, originate mainly in Type la supernovae. If this 
abundance ratio, denoted a/Fe, differs from the solar value, 
then the models may show discrepan cies with the data. A 
remedy for this problem suggested bv lThomas et all (|2003l ) 
is to combine Mg and Fe indices in varying proportions to 
give a composite index that is dependent on metallicity, but 
insensitive to a enhancement. 

In this study we consider the [MgFe]' index defined by 
iThomas et al] (|2003l ). 

[MgFe]' = v/Mgi,(0.72Fe5270 + 0.28Fe5335) (1) 

and the two composite indices of lBruzual fc CharlotI l|2003l ). 

[MgiFe] = 0.6Mgi + 0.4 log(Fe4531 + Fe5015) (2) 

[MgjFe] = 0.6Mg2 + 0.4 log(Fe4531 + Fe5015) (3) 

The latte r two indices were designed to provide a good 
match to iBruzual fc CharlotI (|2003l ) spectra for a broad 
range of galaxy types, but initial studies here showed that 
they do not give as good a fit to the LRG spectra compared 
with our PEGASE models. For this reason we propose two 
additional composite indices, designed to be independent of 
the abundance ratio, whilst also providing observational val- 
ues consistent with those measured for the LRGs using the 
PEGASE model spectra: 

[MgiFe]' = O.SMgi + 0.2 log(Fe5709 + Fe5782) (4) 

[Mg2Fe]' = 0.8Mg2 + 0.2 log(Fe4383 + Fe5406) (5) 

Using data from lThomas et al.l (|2003l ). Figure [1] shows 
that the two indices are, in fact, insensitive to changes in 
Q enhancement^ For both indices, this behaviour is almost 
independent of age or metallicity. In the case of [Mgj^Fe]', 
a/Fe insensitivity holds for models older than ~ 1 Gyr, and 

^ The data used in the plot are available electronically at www- 
astro, physics . ox. ac. uk/ ~ dthomas 
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Figure 1. Absorption indices as a function of a/Fe at fixed over- 
all metallicity. Shown are the new [MgjFe]' (top) and [Mg2Fe]' 
(bottom) indices with their co nstituent Mg and Fe components. 
Data for the plots comes from [Thomas et al. and is rep- 

resentative of_a_J^_G2rold_SSP model with solar metallicity (cf. 
figure 7 of lThomas et al.l ||2003| ')'). 

provides a good metallicity indicator. The [MgjFe]' index is 
less well behaved, with a enhancement becoming apparent 
for models younger than ~ 4 Gyr, especially for low metal- 
licity populations. 

To successfully extract age and metallicity information 
from the SDSS spectra, we choose to simultaneously fit age- 
sensitive and metal-sensitive indices that are both well re- 
produced by our models and that show a weak a/Fe de- 
pendence. After extensive studies fitting a variety of index 
combinations from PEGASE to those from the LRG sam- 
ple, we settle on the following set of indices to be used in 
the final analysis. The indices [MgFe]' and [Mgj^Fe]' are cho- 
sen for their sensitivity to metallicity. It was decided that 
[MgjFe]' should not be used as it is too dependent on a/Fe 
and does not provide an adequate fit to our models. For 
age sensitive indices we choose H/3, H7F and D„(4000). De- 
spite the fact that H7F h as been shown to be sensitive to 
element abundance ratios (jThomas et alj |2003| : iKorn et al.l 
I2OO5I ). including it does not systematically affect our results. 

3.3 Star formation library 

We adopt a Bayesian statistical approach to derive the star 
formation histories, metallicities and masses for each galaxy 



in our sample. The method is similar to the one described 
in K03 and G05, where the data are assumed to be drawn 
from a distribution of models, each described by a parameter 
vector P. A likelihood distribution of a given parameter may 
be obtained by comparing the observational data with a set 
of models that populate the space of all possible P. This set 
of models follows a prior distribution that represents any 
knowledge of various P values in the absence of any data. 

As in K03, we generate a library of 31000 Monte Carlo 
realizations of different star formation histories. Our mod- 
els are characterised by an underlying model where stars 
have been forming according to the exponential law SFR oc 
exp(— 7t) since time tform. Random bursts are then superim- 
posed on this continuous model, and can occur at any time 
after ttorm with equal probability. 

Some of our parameter choices differ from those of K03 
and G05 both because our analysis is restricted to LRGs 
rather than the broad range of galaxy types considered in 
these earlier studies, and to allow a wider possible range in 
the prior. We take tform to be uniformly distributed between 
13.5 and 1.5 Gyr and 7 over the interval to 2. This dif- 
fers from G05, who take 7 between and 1. We also set the 
fraction of models that have experienced a burst in the last 
2 Gyr to 1 per cent, compared with 20 per cent used by K03 
and 10 per cent used by G05. The implications of changing 
these priors is dicussed in Section [4.41 During a burst, stars 
can form at a constant rate for a time tburst, uniformly dis- 
tributed from 3 X 10'' to 3 X 10** yr. The mass of stars formed 
in a burst, Mburst, is characterised by A = Mburst/Mcont, 
where Mcont is the total mass of stars formed by the con- 
tinuous model. A is distributed logarithmically in the range 
0.03 to 4.0. The metallicity of the models is logarithmically 
distributed from 0.02 — 2.5Zq and is permitted to take on 
different values for the continuous model and burst. In most 
cases, however, for simplicity we shall quote the overall light- 
weighted metallicity of the model. 

Multiple versions of each model are then constructed 
by gaussian smoothing the spectra to mimic the stellar 
velocity dispersion of a galaxy. In the range 150 < a < 
300 kms~^, spectra are smoothed at 10 kms~^ intervals, ris- 
ing to 20 kms"^ for 300 < a < 450 kms"^ These intervals 
are chosen so that over 95 per cent of LRG spectra have a 
velocity dispersion measurement that lies within one stan- 
dard deviation of a model value. (The mean velocity disper- 
sion error for the sample is always larger than 20 kms~^.) 
This gives us a total of 23 libraries, each smoothed with a 
different velocity dispersion and each containing 31000 star 
formation histories. 

For each spectrum in these libraries, we compute the 
strengths of the H/3, Rjf, D„(4000), [MgFe]' and [Mg^Fe]' 
indices. The apparent u, g, r, i and z magnitudes are calcu- 
lated using SDSS filters at redshifts between 0.15 and 0.4 in 
intervals of 0.01. Unless otherwise stated, we assume a fiat 
Friedmann-Robertson- Walker cosmology, with Qm = 0.25, 
= 0.75 and Ho = 73 kms'^Mpc'^ 

3.4 Parameter estimation 

The probability density functions (PDFs) of various physi- 
cal parameters are calculated by comparing our set of index 
strengths for a single galaxy to every model spectrum in the 
library. First of all we select the library with the velocity 
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dispersion closest to the observed value. Each model is as- 
signed a weight w = exp(— where is calculated as 
follows: 



X 



E 



o-i 



(6) 



where 7™°'^, 1°^^ and aj are the model index, observed index 
and error values respectively. The sum runs over all five 
indices in our set. The PDF of a particular parameter is then 
given by the distributions of the weights w of all the models 
in the library. We characterize these distributions using the 
mode (the peak of the distribution) and the median. Errors 
are often quoted using the 16 — 84 percentile range of the 
PDF, corresponding to ±1(t for a gaussian distribution. 

In Figure (2] we show some examples of spectra from 
the LRG sample, with the model spectrum that gives the 
best fit based on our set of five indices. In general the best 
fitting spectrum reproduces many of the absoption features 
correctly. However, as we only fit absorption indices, the 
overall shape of the continuum is less well reproduced, espe- 
cially towards the red end of the spectrum. The strength of 
this method is that we can derive the star formation history 
independently of the shape of the continuum and the extinc- 
tion model. We only need an estimate of the dust extinction 
when computing mass estimates. 

To characterize the age of each galaxy in our library, we 
calculate the bolometric light-weighted age 



tL = 



J* SFR{t - T)L{T)TdT 

jl SFR(t - T)L{T)dr ' 



(7) 



where L{t) is the bolometric luminosity of the system 
evolved over a time r. Similarly, we can define a light- 
weighted metallicity 



Zl = 



/J SFR{t ~ T)L{T)Z{t - r)dr 
/J SFR{t - r)L{T)dr 



(8) 



where Z{t~T) is the metallicity of the stars at the time t — r. 
The metallicity may also be written as [Z/H] = log{Z/ZQ), 
where Zq = 0.02. In this study we will refer to these quan- 
tities as the age and metallicity of a model, respectively. 

To estimate the stellar mass of a galaxy, we first need 
to calculate the effect of dust on the spectrum. Here we 
model the dust as a simple slab of material between us and 
the galaxy. Any light, ^(A), from the galaxy is attenuated 
according to the following wavelength-dependent law, 

Sob.(A) = Semit(A) X lO-"-*^-*^ (9) 

where A(A) is given by the extinction curve o f ICardeUi et al.l 
(|l989[ ). This empirical model takes the following two- 
component form: 

where a(A) and 6(A) are pre-determined polynomials. We 
adopt the standard Milky Way value of Rv =3.1 when fit- 
ting dust models to the data. This gives essentially a single- 
component dust model, characterised by Ay, the attenua- 
tion in the V-band (centred on 5500 A). We infer the value 
of Av using the difference between the fibre r — i colour and 
the r — i colour of each model at a redshift matching the 
data, 
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Figure 3. The distribution of dust attenuation values in the z 
band for the 4391 LRGs in our sample. 
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Using equation 1101 we then calculate the extinction- 
corrected, photometric z-band luminosity of the entire 
galaxy. The stellar mass for a given model can be estimated 
as follows. 



X 



(-) 



(12) 



where (M,/L2)niodci is the redshifted mass-to-light ratio in 
the 2-band. This assumes that the mass-to-light ratio within 
the fibre is the same as for the overall galaxy. G05 find no 
significant trend of age or metallicity with the fraction of 
total light entering the fibre for massive early-type galax- 
ies, so this seems a fairly safe assumption. The mass values 
for each model are then weighted as usual to give a PDF 
for the stellar mass. As in G05, we allow attenuations down 
to Az = —0.1 to account for errors in measurements and 
discrepancies within the models. Figure |3] shows the distri- 
bution of median Az values for our sample of SDSS LRGs. 
The typical attenuation in the z band is quite small, with 
a median value of 2.6 x 10"'^. This is consistent with the 
picture of LRGs as early-type galaxies containing old stellar 
populations with little dust and gas. 



4 RESULTS 

We now apply the above analysis to the SDSS sample of 
4391 Luminous Red Galaxies as described in Section [2l 



4.1 Index fitting 

Figure|4]shows the distribution of indices in the LRG sample 
when fitted with the final set of indices , using the method 
described in Section [3.31 (cf.. figure 18 of iBruzual fc CharlotI 
(2003')). While this plot shows that many indices agree very 
weU with models (e.g. [MgFe]', [Mg^Fe]', Ca4227, Fe5782, 
H/3, H7F, Dn(4000) ) others show comp lete failure (e.g. 
C24668, NaD, US a). iThomas et all (|2003l l and iKorn etahl 
(2005) investigate the effect of varying individual ele ment 
abundances on the Lick indices. iThomas et all l|2003l ) find 
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Figure 2. Examples of LRG spectra from the SDSS are shown in green, after dereddening for Milky Way extinction. The PEGASE 
spectra from the Monte Carlo library with the best-fitting set of the selected five absorption indices (see text) are overplotted in red. 
Top-Left: The best fit shows a young population with a light-weighted age of 3.5 Gyr and metallicity Z = 0.03. Top-Right: An example 
of an old population that formed 9.9 Gyr ago, which experienced a significant burst 1.7 Gyr ago. The best-fit metallicity is Z = 0.02. 
Bottom-Left: An example of a spectrum that our method finds a good fit for despite a low SNR (6.3). The galaxy is 5.7 Gyr old with a 
metallicity of Z = 0.03. Bottom-Right: An example of a spectrum for which our library fails to find a good fit. The best fitting set of 
indices gives = 19.7 for five degrees of freedom. The previous three spectra have = 4.7, 5.2, 2.0 respectively. The model spectra 
have been smoothed slightly for presentation. 



that model predictions for C2 4668 poorly reproduce glob- 
ular clust er data, eve n when Mgi, CNi and CN2 are well- 
modelled. iKorn et al. find this index is most sensitive 
to the carbon abundance. They also find that the NaD in- 
dex is sensitive to both sodium and metallicity and H^a is 
affected by iron abundance and metallicity. The narrower 
passband of H5f reduces the sensitivity to iron. The blue 
indices CNi, CN2, Ca4227, G4300, Ca4455, Fe4531 and 
C24668, along with Ti 02, are known to suffer the worst 
modelling uncertainties (|Korn et al.ll2005l ). Many of the re- 
maining discrepancies are likely caused by overabundances 
of individual elements that are not accounted for by the 
PEGASE model. For instance, the iron index Fe5015 is sensi- 
tive to Ti and Mg in additio n to iron, and to a enh ancement 
as well l|Thomas et al.ll2003l ). lThomas et all (120031 ) also find 
that the Mg;, index is sensitive to a enhancement, more so 
than are either the Mgi or Mg2 indices, for which the agree- 
ment is better, though the spread is still large. The under- 
prediction of some indices, li ke CNi and CN2, is expected 
for a-enhanced abundances l|Thomas et al.l [20031 ). We dis- 
cuss the indices further in Section [S] 



4.2 LRG properties 

Here we present the results obtained for the light-weighted 
age, metallicity and stellar mass of the 4391 Luminous Red 
Galaxies in our sample. Figure [5] shows the distributions of 
the median-likelihood estimates of these parameters (left- 
hand side) and their errors (right-hand side) in the top, 
middle and bottom panels respectively. The age distribu- 
tion shows a peak at around 6 Gyr, suggesting the LRGs 
are comprised of mainly old stellar populations. Typical un- 
certainties on the age are (Tt ~ 2 Gyr, with no galaxies show- 
ing errors larger the 3 Gyr. The distributions also indicate 
an average super-solar metallicity of [Z/H] ~ 0.2, with a 
long tail down to [Z/H] ~ —0.6. Errors in the metallicity 
measured in individual galaxies are around (J[z/h] ~ 0.15, 
extending up to 0.3 for some galaxies. Similar distributions 
are found for the age and metallicity of early type galaxies by 
iTrager et al.l l|2000bl ). Finally, the mass shows a narrow dis- 
tribution that peaks at log{M, /Mq) ~ 11.5. It is fairly well 



constrained, with a maximum error of ai, 



og(Af./A/Q) 



0.1. 



The high mass of the sample is expected from the LRG high 



8 Tom Barber, Avery Meiksin and Tara Murphy 



T — n — I — r 

[MgFe]' 



1 — I — I — I — r 

[Mg,Fe] 



1 — I — I — I — r 

[MggFe] 



T — I — I — I — r 

[Mg^Fe]' 



1 — I — I — I — r 

[MggFe]' 



0.6 
0.4 
0.2 
0.0 
0.6 
0.4 
0.2 
0.0 
0.6 
0.4 
0.2 
0.0 

■n 

I 0.6 

o' 

o 0.4 

(O 

§ 0.2 

CD 
V) 

0.0 
0.6 
0.4 
0.2 
0.0 
0.6 
0.4 
0.2 
0.0 
0.6 
0.4 
0.2 
0.0 




HTa+H8a 







Ca4227 



G4300 








Fe4383 



Ca4455 



Fe4531 



Cp4668 






Fe5015 



Mgg 



Fe5270 







Fe5335 



Fe5406 



Fe5709 



Fe5782 



NaD 








TiO, 



TiO, 







D(4000) 



HKratio 




J. 




-4.0 -2.0 0.0 2.0 4.0 -4.0 -2.0 0.0 2.0 4.0 -4.0 -2.0 0.0 2.0 4.0 -4.0 -2.0 0.0 2.0 4.0 -4.0 -2.0 0.0 2.0 4.0 

(l'"-l*=)/a| 



Figure 4. Simultaneous fits of the H^, Hjf, D„(4000), [MgFe]' and [MgjFe]' indices (highUghted) of the 4391 galaxies from the SDSS 
LRG sample with measurable indices. Each panel shows the distribution of the best fitting index minus the observed one, 7°*"*, 
divided by the index error crj. The dashed line in each panel indicates a gaussian with unit standard deviation. 
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Figure 5. Histograms of tlie median values of age (top-left), 
metallicity (middle-left) and stellar mass (bottom-left) for 4391 
Luminous Red Galaxies from the SDSS-DR4. The y-axis units 
give the fraction of galaxies in a particular bin. The right-hand 
panels show the respective error distributions, calculated as one- 
half of the 16-84 percentile range. 



luminosity cuts. The narrow distribution suggests that we 
have indeed selected a homogeneous population of galaxies. 



4.3 Age, mass and metallicity relationships 

We are also interested in exploring relationships between 
the mass, age and metallicity of the galaxies in our sample. 
Rather than simply plot the median of two separate distri- 
butions, we calculate the joint distribution of each pair of 
parameters. This is done by adding the 2D normalised prob- 
ability distribution of all galaxies in our sample, giving each 
galaxy equal weight. 

In this study, we use the same cosmology as that of 
iDe Lucia et al.l (|2006| ). who use merger trees from the Mil- 
lenium Simulation, grafted onto semi-analytic models of 
galaxy formation, to track the properties of elliptical galax- 
ies. As LRGs are early-type galaxies, we can directly test the 
predictions of these simulations at the high-mass end using 
the LRG sample. 

In the top- left panel of Figure (6] we show the 2D prob- 
ability distribution of age against mass. In this plot, the age 
has been estimated for each galaxy at a redshift of z = 0, 
by adding the lookback time to the light-weighted age. Also 
shown is the distribution of metallicity against stellar mass 



in the m iddle- left. These two plots are comparable to fig- 
ure 6 of iDe Lucia et al.l l|2006l ) , whose predictions for the 
same variables are shown as the solid and dashed lines, rep- 
resenting the median and upper and lower quartiles respec- 
tively. 

The plots do not show any strong trends with stellar 
mass, and reflect the distributions seen in Figure [5] The 
age appears to be slightly correlated with mass, and sug- 
gests that the less massive galaxies are made up of younger 
stellar populations within the restricted mass range of the 
LRGs. The metallicity appears to be fairly uniform over 
the stellar mass range studied, with at most a slight anti- 
correlation with mass. The slight bimodality in the distri- 
bution is discussed further in Section 14.51 The ages pre- 
dicted bv lDe Lucia et al.1 (|2006l ) exceed our estimates by 1- 
3 Gyr. Their metallicity predictions agree closely with ours, 
although their predicted values are lower by about 0.1 dex. 

An issue that invariably arises when deriving galaxy 
properties is the age-metallicity degeneracy, which oc- 
curs due to competition between age and metallicity on 
absorption- line strengths, so that old, metal-poor galaxies 
can be mistaken for young, metal-rich ones (or vice versa). 
Figure |6l bottom-left panel, gives the 2D likelihood distri- 
bution of metallicity against light-weighted age (adjusted to 
z — 0), with the degeneracy manifesting itself as an apparent 
anti-correlation between the two parameters. The bulk of the 
distribution is centred on ages of 6-7 Gyr and [Z/H]~0.2, 
reflecting the ID distributions in Figure [5] To quantify this 
degeneracy we calculate the median metallicity in each of the 
age bins of Figure |S] and fit a straight line through the result- 
ing points. This gives a slope of A log Z/ A log t — —1.05 dex 

Previous studies (ICaldwell et al.l l2003l : iKelson et al.l 

have suggested possible correlations between the age 
and me tallicity of ear l y-typ e galaxies with velocity disper- 
sion, a. iKelson et al.l (|2006l ) go so far as to claim that the 
properties of early-type galaxies in clusters can be com- 
pletely characterised as a one-parameter family, given by 
the velocity dispersion of the galaxy. The Faber-Jackson re- 
lation also motivates investigation into the dependence of 
mass on velocity dispersion. The right-hand column of Fig- 
ure |S] shows how the three parameters, age, metallicity and 
mass, depend on velocity dispersion in our sample. The age 
seems to show a small correlation with velocity dispersion: 
the median yields a slope of Aage/Aa ~ 4 Myr/kms~^. The 
metallicity shows an anti-correlation with A[Z/ H]/ Aa ~ 
— 1 X 10"'^ dex/kms~^. The mass shows a slight positive cor- 
relation of Alog(M/MQ)/Acr ~ 9 X 10"'' dex/kms"^ The 
modes of the distribution, however, are well concentrated, 
centred at about 240 km s"'^. 



4.4 Star formation history of LRGs 

In addition to deriving the rest-frame parameters of these 
galaxies, we would like to constrain their star formation his- 
tories by tracking the overall star formation rate as a func- 
tion of cosmic time. To achieve this, we compute the look- 
back time for each galaxy and extrapolate the past and fu- 
ture star formation histories for each model, weighted in the 
usual way. This gives us a PDF for star formation rate, in 
bins of lookback time, which can be summed over the entire 
sample to give an overall likelihood distribution. It should be 
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Figure 6. Two-dimensional probability distributions for various pairs of galaxy properties. The plots were obtained by adding the 
normalized 2D distributions for all 4391 galaxies in the SDSS-DR4 sample of LRGs. The left hand-column shows light-weighted age (at 
2 = 0) and metallicity against stellar mass, and metallicity against age, from top to bottom respectively. The right-hand column shows 
age, metallicity and mass against vel o city d ispersion. The solid line in the upper two plots of the left-hand column shows the median 
value as predicted by iDe Lucia et al] ||2006| ). using semi-empirical models. The dashed lines represent the upper and lower quartiles of 
their distribution. The solid lines in the right-hand column plots give the median parameter value in each velocity dispersion bin. 
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Figure 7. The median star formation rate of LRGs in our sample, 
binned by lookback time. (See text for details.) The continuous 
lines are the semi-analytic model predictions of iDe Lucia et al] 

noted that this is only done for the continuous component 
of the models, and does not include star formation from the 
burst. We also note that the earlier studies of K03 and G05 
used Monte Carlo libraries generated with < 7 < 1 for 
the continuous component, -while here -we choose < 7 < 2. 
We discuss belo-w the impact of changing this prior on the 
derived star formation history. 

In Figure [T] we present the median values of star for- 
mation rate per stellar mass, binned by lookback time. The 
error bars show the 16th and 84th percentile limits of the 
distribut ions in each bin . Also shown is the formation rate 
from the lDe Lucia et al] (|2006l ) models (from their figure 1) 
with stellar masses of 10^^ M© and 10^^ A/©, for galaxies 
in clusters and in the field. For reference, the LRG sample 
used here is selected in the redshift range 0.15 < z < 0.4, 
corresponding to lookback times of 1.8 Gyr < tit < 4.2 Gyr. 

Between lookback times of 3 and 12 Gyr, the LRG 
results and the simulation show good agreement, with all 
points lying within the 68-percentile range of the predicted 
value. For shorter lookback times, the data give a lower star 
formation rate than expected when compared with simu- 
lated field galaxies. In this regime of lookback time, we are 
extrapolating the future star formation of a galaxy using 
the passive exponential model and take no account of any 
possible interactions. This may account for the discrepancy, 
as additional mergers may boost the star formation rate. In- 
terestingly, there is much better agreement at more recent 
times when considering galaxies in clusters. This is consis- 
tent with the inference that m ost LRGs reside in clusters 
(|Kormendv fc Diorg:ovskilll989l ). At lookback times exceed- 
ing 12 Gyr, the model curves drop off sharply, while the data 
still reveal some residual star formation. This is because we 
only show star formation rates for galaxies with stellar pop- 
ulations existing at these early times. Accordingly, the mea- 
sured star formation rate at 12 Gyr accounts for only a small 
fraction of galaxies that will be created in the Universe. 

When considering the overall formation mechanism of 
elliptical galaxies, it is important to distinguish between the 
formation times of the composite stars and the assembly 
time of the constituent progenitors. The formation time of 
a system can be defined as the time when a certain fraction 
of stellar mass has been formed, compared with the mass at 




Redshift Lookback Time / Gyr 

Figure 8. The redshift (left-hand column) and lookback time 
(right-hand column) at the formation time and when 50 per cent 
and 80 per cent of the stellar mass is in place, compared with 
the mass at 2 = 0. (Top, middle and bottom panels respectively.) 
The solid line indicates that distribution for a library in which 
0.0 < 7 < 2.0, while the dotted line corresponds to 0.0 < 7 < 1.0. 
The lower two distributio ns (50 per cent and 80 per cent) are 
comparable to Fig ure 4 of iDe Lucia et al.l 1I2OO6I '). 



z — 0. Here we calculate the cosmic lookback times corre- 
sponding to the redshift at which star formation first begins, 
and when 50 per cent and 80 per cent of the stellar mass has 
been formed, denoted Tform, Tf50 and Tf80 respectively. 
These parameters give us not only a handle on when the 
stars in the galaxies formed, but the time-scale over which 
the stars which reside in LRGs today were forming. 

The solid histograms in Figure [S] show the distributions 
of these formation times and their corresponding redshifts. 
Table [1] summarizes these results, giving the lower quartile, 
median and upper quartile values. The results show that 
the stellar component of LRGs typically first began forming 
around 8-10 Gyr ago (2 « 1.1 — 1.9), with 80 per cent of the 
stellar mass in place 6-8 Gyr ago {z ~ 0.7 — 1.1). 

These values are directly compa r able w ith those pro- 
vided in table 1 of iDe Lucia et al.l (|2006l ). whose semi- 
analytic models predict Tf50 ~ 11 Gyr and Tf80 ~ 9 - 
10 Gyr. Our results for < 7 < 2 give ages about 2-3 Gyr 
younger than the simulation predictions. 

In Figure [5] (top), we compare the distribution of the 
median 7-values when < 7 < 2 is allowed and when the 
library is restricted to models for which < 7 < 1. The 
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25% 


50% 


75% 


< 7 < 2 


Tform 
Tf50 
Tf80 


8.15 
7.52 
6.63 


9.27 
8.58 
7.63 


10.59 
9.97 

9.12 


< 7 < 1 


Tform 
Tf50 
Tf80 


10.38 
9.37 
8.03 


11.08 
10.17 
8.90 


12.14 
11.22 
9.74 


De Lucia et al. ('2006) 
(10.8 < log(iW/M0) < 11.5) 


Tf50 
Tf80 


9.78 
7.56 


10.56 
9.20 


11.22 
10.31 


De Lucia et al. (2006) 
(11.5 < log(M/MQ) < 12.2) 


Tf50 
Tf80 


11.01 
9.50 


11.41 
10.31 


11.76 
10.79 



Table 1. The lower quartile (25 per cent), median (50 per cent) 
and upper quartile (75 per cent) values summarizing the lookback 
times (in Gyr) since formation of LRGs in our sample. Values are 
shown for the full Monte Carlo library, for which < 7 < 2, 
and also the results when the library is restricted to models with 
< 7 < 1. For referenc e, the corresponding values from table 1 
of lDe Lucia et al.l l|2006l ) are reproduced here for comparison. 

diflterence in predicted values is striking, with median val- 
ues of 1.28 and 0.73 for the two distributions, respectively. 
The impact that this has on the predicted star formation 
history is shown in Figure |9] (bottom). This plot gives the 
2D probability distribution of 7 as a function of the time 
since formation of the galaxy, produced in the same way as 
Figure (6] We see here that the two parameters appear to 
be inversely proportional to one another, with older objects 
having the most prolonged periods of star formation. Gen- 
erally, only objects older than 7 Gyr are expected to have 
7 < 1. The overall effect of restricting the prior on 7 to 
< 7 < 1 is therefore to increase the derived ages of the 
galaxies. 

This is illustrated by the dashed histograms of Fig- 
ure [8] and the corresponding values in Table [1] Here we see 
that for the more restricted range in 7, the predicted me- 
dian formation time is increased by 1.8, 1.6 and 1.3 Gyr 
for Tform, Tf50 and Tf80 respectively. The inter-quartile 
range for the formation time distributions are reduced from 
2.4 Gyr (0 < 7 < 2) to 1.7 Gyr (0 < 7 < 1). This is 
understandable as by reducing the range of 7, we are also 
reducing the number of possible star formation histories for 
a given formation time. Using the more restricted range of 
7 brings our predictio ns in closer agreement with those of 
iDe Lucia et al.l l|2006l) . We note that restricting the prior 
value of 7 also increases the median light-weighted age by 
1-2 Gyr, but appears to have little effect on the derived 
metallicity and mass distributions. 

In addition to tracking the formation times of galaxies, 
it is important to discover whether galaxies have undergone 
any recent bursts of star formation activity. For each galaxy, 
we calculate the probability Pburat of a burst by adding the 
weights of every model that has experienced a burst in the 
last n Gyr, and normalising by the total weight of all mod- 
els in the library. The distribution of burst probabilities for 
n = 2, 4, 6 is shown in Figure fTOl (bottom panel). The proba- 
bilities are generally small. Typically there is only a 0.25 per 
cent chance of a burst occuring within the last 2 Gyr for an 
individual galaxy. This rises to 1.6 per cent for bursts within 
the last 4 Gyr and 3.2 per cent within 6 Gyr. We note that 
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Figure 9. The top panel shows the distribution of 7 for the 
LRG sample when the prior range is 0.0 < 7 < 2.0 (solid) and 
0.0 < 7 < 1.0 (dashed). The lower panel shows the 2D PDF for 
7 as a function of the lookback time corresponding to the time of 
formation. 

the burst probabilities may be sensitive to the prior. For 
reference, we set up the library so that 1 per cent of galax- 
ies experience a burst within the last 2 Gyr. We would also 
like to calculate the magnitude of these bursts compared 
with the underlying continuous model. To do this we con- 
sider only those models in which bursts have occured within 
the lifetime of the galaxy and calculate the median value of 
A — Mburst /A^cont by marginalising in the usual way. The 
top panel of Figure [TD] gives a typical value of A = 0.4, im- 
plying that around 29 per cent of the total mass of stars in 
a galaxy that had a burst in its history would have come 
from the burst. 

4.5 Metallicity of LRGs 

In this section, we discuss implications of the measured 
metallicity distribution of the LRGs in terms of their 
chemical evolution. The simplest possible model for the 
chemical evolution of a galaxy is the closed-box model of 
iTalbot fc Arnettl(|l97ll ). We also consider a standard modi- 
fication allowing for continued accretio n of gas by the galaxy. 
Detail s of both models are discussed in iBinnev fc Merrifieldl 
(199^. 

In the closed box model, it is assumed that initially the 
material is entirely gaseous, possibly enriched to a metallic- 
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Figure 10. The top panel shows the distribution of A, the ratio of 
mass formed in a burst compared with that formed by the expo- 
nentially decaying model. The lower panel shows the distribution 
of burst probabilities. The solid, dashed and dotted lines show the 
distribution of Pbursti the probability that a galaxy experienced 
a burst, in the last 2, 4 and 6 Gyr, respectively. 



ity level Zi. As the name suggests, no material is allowed to 
enter or leave the region, so there is no further appreciable 
interaction of the galaxy with its surroundings once it has 
formed. As the galaxy evolves, stars are formed, consuming 
interstellar gas in the process and polluting the remaining 
gas with heavy elements corresponding to a yield p of mass in 
metals per mass of stars formed. The yield is assumed con- 
stant. It is further assumed that turbulence keeps the gas 
homogeneous so that the instantaneous recycling approxi- 
mation may be made, which neglects the delay between the 
formation of a population of stars and the ejection of heavy 
elements into the remaining gas by massive stars in the pop- 
ulation. 

Suppose the interstellar gas has a mass Aig and con- 
tains a mass of heavy elements Mh- The metallicity of the 
gas is then 

In this model, the metallicity may be shown to evolve ac- 
cording to 



Zit) 



■ pin 



(14) 



MM 

The mass of stars in the box with metallicity between Z,, 



and Z{t) is then 

Msi<Z) = Ms{t)=Mg{0)-Mg{t) 

Z - Zi 



MgiO) 



exp 



P 



and the metallicity distribution of the stars is 



dMs{< Z) _ Mg{0) 



dZ 



P 



exp 



(15) 



(16) 



for Z > Zi, and vanishes for 
of the stars is then 

rZ(t) 



Z < Zi. The mean metallicity 



,dMs 
' dZ 

.T In X ^ 



dZ 



/ , .T In X \ 

Ms{t)/Mg{0) 



(17) 



where x = 1 - Ms{t)/Mg{0) = exp[-(Z - Zi)/p]. The 
limiting metallicity values of Z = Zi and Z ^ oo corre- 
spond respectively to a; = 1 and a; = 0, and to Z — Zi and 
Z — Zi+p. The mean metallicity is thus restricted to being 
smaller than the sum of the initial metallicity and the yield. 
Inspection of Figure [TT] suggests the yield must then be as 
high as p ~ 0.05. 

In this scenario, the distribution of mean metallicity 
is a reflection of the distribution of initial metallicities Zi 
of the galaxies and their ages, with younger galaxies ex- 
pected to have lower metallicities on average. However, this 
picture is not supported by the age-metallicity plot of Fig- 
ure [6] which, while showing a large amount of scatter for 
super-solar metallicity values, also shows a general anti- 
correlation, particularly at low metallicity values. This effect 
is likely in part due to a residual age-metallicity degeneracy, 
although the degree to which the degeneracy masks a true 
underlying trend is unclear. 

We next consider a modification of the closed box model 
in which star formation is abruptly halted after an enrich- 
ment phase. The model is motivated by the proposed solu- 
tion to the anti-hierarchical behaviour of early type galax- 
ies invoking a sudden quenching of star formation in the 
sub-units which eventually coalesce into the galaxies we see 
today. In this model, we allow for stochastic star-formation 
quenching (SSQ), modelled as a Poisson process at a uni- 
form rate A. The probability that star formation is suddenly 
quenched in any given galaxy between the times t and t + dt 
is then Xdt. The metallicity level any galaxy achieves is then 
given by the time during which the gas was enriched before 
star formation ceased, as well as by any initial metallicity. 
The metallicity distribution is given by 

1 dAf -xt.dt 
—z ^ = e A — ^. (18) 

Mot dz dz ^ ' 

The time derivative dZ jdt requires specifying a star forma- 
tion rate ?/'(t). The simplest approximation is a constant rate 
■!/j(t) = ■i/'o- For reasons explained below, we instead consider 
the more general form 

-0 



(19) 



where is a constant and /3 < 1. A constant rate 
corresponds to the special case /3 = 0. For /3 > 0, the rate 
corresponds to a burst at t = 0, gradually decaying after- 
wards. The restriction /9 < 1 ensures a finite mass of stars 
is formed. 
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Expressing dt/dZ 
ip{t) — dMa/dt gives 



x{t) = 1 - 

and 

p dM 
A/tot dZ 



V'/3 



[dt / dx) / [dZ / dx) and noting 



(20) 



1-/3 



exp [-a(l-a;)i/(i-«] 



2-/3 

(I-X)T^ 

1 — x + In x 



(21) 



where a = A[A4g(0)(l -/3)/V'/3]^^*^'''' is the ratio of the gas 
consumption time-scale to quenching time-scale. Typically 
this ratio is expected to be large, otherwise most of the mass 
will have been enriched to its maximum mean value through 
generations of star formation before star formation ceases. 

In the limit of small mean metallicity values, Z 
Zi, the distribution function (p / J\ftot)dJ\f / dZ 2y^(1 — 
j,^/3/(i-/3) {-Qj. 2; ^ 1. For the special case P — 0, the limit- 
ing value is 2a. This limiting value will in fact apply to any 
analytic form for ip{t) that is non- vanishing at t = 0, in- 
cluding an exponentially decaying star formation rate. Such 
asymptotic behaviour, however, is contrary to the diminish- 
ing distribution at low metallicity values in Figure 1111 For 
< /3 < 1, the distribution function vanishes for Z Zi. 
This is the reason for having introduced this particular 
'cuspy' form for the decaying burst. 

We next modify the simple closed box model by as- 
suming that metal free gas may be accreted by the galaxy 
from its surroundings. The change in total mass is SMt = 
5A4s + SMg 7^ 0. In this case, the metallicity evolves ac- 
cording to 



Z{t) = p 



1-11-^ 
P 



(22) 



where x — 1 — Ms{t) /Mg{Q), as for the closed box case. The 
relation between x and t is again given by Eq. 1201 As there 
is no upper limit to the amount of gas accreted, there is no 
upper limit to the amount of stars which may form. As a 
consequence, x may become arbitrarily small. In the limit 
X — > —00, the metallicity takes on the value of the yield 
{Z ^ p). The maximum metallicity achievable is the yield, 
provided Zi < p. (For Zi > p, the metallicity will decrease 
with time towards p.) 

The mean metallicity evolves according to 



Z{t) = p - 



Z{t) - Zi 
ln£^ ' 



(23) 



with the limiting value Z —> p as x ^ —oo. In analogy to 
the closed box case, the metallicity distribution in the SSQ 
model is given by 

dj\r 



_P ^ 

A/tot dZ 



1-/3 



exp [— Q!(l — X 
(1 - 2;)T^ 



1/(1-/3) 



+ 1 



Z-Zj 



+ Zr-p 



(24) 



Best fits to the measured distribution are shown in Fig- 
ure [TT] (upper panel) for both the closed box (solid line) 
and accreting box (dashed line) SSQ models. For the closed 
box SSQ model, the parameters for the fit are Zi = 0.0040, 
Q = 66, /3 = 0.59 and p — 0.26. For the accreting box 
SSQ model, the parameters for the fit are Zi = 0.0040, 



a = 2.3, (3 = 0.50 and p = 0.10. Both models well describe 
the distributions. Error bars are not quoted because neither 
fit has an acceptable value of (reduced X'^ 5), due in 
large part to non-Poisson fiuctuations from bin to bin. A dip 
near solar metallicity hints at a bimodal distribution. With- 
out extensive Monte Carlo tests of the population synthesis 
model- fitting procedure, it is difficult to assess the statis- 
tical significance of the fluctuations, and of any evidence 
for bimodality. Doing so is beyond the scope of this paper, 
but may be worth considering in any future estimates of the 
galaxy metallicity distribution. The spread in parameter val- 
ues between the two models provides some indication of the 
errors. We note that the excess number density of systems 
predicted by the models for Z ~ 0.05 may be an artefact of 
the stellar library, which is restricted to Z ^ 0.05. 

Values for a exceeding unity are expected, and con- 
firm the assumption that star formation is quenched on 
a time-scale short compared with the time-scale for com- 
plete conversion of the initial gas mass into stars. Adopting 
two e-foldings of star formation in the population synthe- 
sis models for the conversion time, the distribution of 7 in 
Figure [9] suggests a typical time-scale for depleting the ini- 
tial gas content of 2/7 ~ 1.5 x 10® yr. The corresponding 
quenching time-scale A^^ for the closed box SSQ model is 
then 2/70: « 2 X 10^ yr, comparable to the Salpeter time 
(isaip ~ eaTc/A-nGm-p) for a black hole energy conversion 
efficiency e of 5 per cent, and to the lower limi t for the 
lifetime of a luminous QSO (|Yu fc Tremaindl2002h . The ef- 
fective yield value for the model is quite hig h, comparable 
to that expected f rom massive stars alone l|Maedeij Il992l : 
iHirschi et allbOOSl ), and so is not expected for a normal 
initial mass function. Such an IMF would result almost en- 
tirely in stars sufficiently massive to produce Type II super- 
novae. A metallicity value corresponding to Type II super- 
novae ejecta is consistent with the evidence for a-enhanced 
abundances. An IMF dominated by massive stars should not 
necessarily be precluded. An initial mass function heavily 
weighted towards massive stars {dN/dm ~ I/jti) has been 
invoked in recent semi-analytic galaxy formation models to 
reproduce sub-millimetre galaxy counts and the high metal- 
licities observed in the intracluster gas of rich galaxy clusters 
(Lacey et al. 2006, in preparation). Luminous Red Galaxies 
are believed to reside predominantly in the cores of just such 
rich clusters. 

For the accreting box SSQ model, the smaller value for 
Q corresponds to a star-formation quenching time-scale of 
~ 6 X 10* yr, com parable to the upp er limit for the lifetime of 
a luminous QSO (|Yu fc Tremainel 20021. The effective yield 
value for the model is still high, about twice that expected 
for the Chabrier (2003) initial mass function for early stars, 
for which a little under half of the stellar mass is contained 
in stars massive enough to form Type II supernovae, and 
assuming an approximate mean ejected mass fraction of 10 
per cent. It suggests again an IMF tilted toward more mas- 
sive stars, but not as strongly as required by the closed box 
model. 

The sudden star-formation quenching scenario may also 
account in part for the inverse correlation between age and 
metallicity. For a given formation time, systems with an ex- 
tended amount of star formation before quenching will pro- 
duce stars with increasing metallicity and a younger average 
population compared with systems in which star formation 
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Figure 11. The top panel shows the metaUicity distribution func- 
tion for LRGs in our sample (histogram). The solid line is a fit of 
a closed box SSQ model, while the dashed line is for an accreting 
box SSQ model (see text). The bottom panel shows the relation 
between age and metaUicity predicted by the accreting box SSQ 
model for characteristic quenching time-scales of 0.65 Gyr (solid 
lines) and 1.3 Gyr (dashed lines) for formation times of 9 Gyr 
and 12 Gyr. 
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was quenched early. Because of the short characteristic decay 
time for star formation, however, the ages of the stars will 
not differ greatly (typically ~ 27"^ ~ 1.5 Gyr), so that the 
inverse relation must also result from a range in formation 
times. For slower star formation, however, a broader range 
in ages is possible. The relation between the lookback time 
to the time at which star formation was quenched and the 
mean metaUicity of the stars at that time is shown in Fig- 
ure lllK lower panel) for the accreting box SSQ model. Values 
of = 0.65 Gyr (solid lines) and — 1.3 Gyr (dashed 
lines) are shown, with formation times of 9 Gyr and 12 Gyr. 
These curves span the range of metaUicity vs age in Fig- 
ure [6] The rapid ascent at low (sub-solar) metallicities, and 
slower metaUicity evolution at higher (super-solar) metal- 
licities, predicted by the model suggests that most galaxies 
will be concentrated in the region of super-solar metallici- 
ties, consistent with the 2D age-metallicity distribution in 
Figure [G] 



5 DISCUSSION & CONCLUSIONS 

In this section, we compare our results with those in previous 
works, and summarise the main conclusions of our study. 
The groups at MPA/JHU have made publically available 
large catalogues of the results of model fits to SDSS data. 
These include the mass catalogues of K03 (updated to Data 
Release 4) and age, mass and metaUicity values of G05 (Data 
Release 20. 

To make our results comparable with values from K03 
and G05, the library is restricted to models with < 7 < 1, 
to more closely match our prior with theirs. We then re- 
compute masses by fitting only the D„(4000) and H5a 
indices to allow direct comparison with K03. Similarly, 
we calculate masses, ages and metallicities by fitting the 
[MgFe]', [Mg2Fe], H/3, YL-^a+^Sa and D„(4000) indices, as 
used in G05. We also match cosmologies, adopting a fiat 
Friedmann-Robertson- Walker universe with Q.m ~ 0.3 and 
Ho = 70 kms-^Mpc~\ 

The left-hand column of Figure [12] shows a compari- 
son between our values and those from the MPA/JHU cat- 
alogue. All derived parameter values show generally good 
consistency over the ranges studied. The right hand-column 
of Figure [12] shows the distribution of differences between 
our values and the MPA/JHU values. The solid histogram 
uses the modified prior and K03 or G05 indices as above. 
The dashed histogram shows the same residuals using the 
full library (0 < 7 < 2) and the original set of indices 
([MgFe]',[MgiFe]', H/3, U'yp and D„(4000)). 

The stellar mass distribution we derive is narrowly 
peaked around 3 x lO^'^ Mq. A comparison between our de- 
rived stellar masses and those of K03 shows overall consis- 
tency, with a scatter of ~ 0.16 dex. Our stellar mass values 
are in general larger than theirs by ~ 0.09 dex when using 
their indices, and ~ 0.07 dex using ours. There is a smaller 
offset when we compare our masses to those derived by G05: 
~ 0.02 dex using their set of indices, and ~ —0.04 dex with 
ours. 

The offsets are likely to be due to a combination of fac- 
tors. Principal among these is the difference in the adopted 
initial mass functi on (IMF). Wh ile we use the Salpeter 
I MF, K03 u se the iKroupal (I2001I) IMF and G05 us e that 
of iChabrieil (|2003| ) for disc stars, feruzual fc" Chariot! (I2003D 
show that, for a simple stellar population with solar metal- 
licity, the Salpeter IMF results in an increased stellar mass- 
to- light ratio of about 0.2 dex over the Kroupa and Chabrier 
IMFs, at a given evolutionary time. Our stellar mass values 
agree with the results of K03 and G05 more closely than 
this, which on the surface is surprising. The reason for the 
good agreement appears in fact to be in part fortuitous. The 
effect of broadening the range of 7 to < 7 < 2 is to find 
ages that are, on average, ~ 1.1 Gyr smaller than those of 
K03 and G05. If we limit the prior to < 7 < 1, this differ- 
ence reduces to ~ 0.4 Gyr. The scatte r in ages is around 
^ 1 G yr. According to the results of iBruzual fc Chariot! 
(l2003h . the larger ages found by K03 and G05 correspond 
to an increase in the mass-to-light ratio of close to 0.1 dex 
over the values corresponding to our best fits. This differ- 
ence brings our stellar mass estimates and those of K03 and 



^ Data available electronically 
garching.mpg.de/SDSS 
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Figure 12. Comparison of results from this paper with previous work. The left-hand column shows our stellar mass values against those 
of K03 and G05 (top two plots), and our age and metallicity against the values of G05 (bottom two plots). The values quoted have been 
calculated by fitting the similar indices and priors as K03 and G05 respectively. The right hand column shows the residual between our 
values and the corresponding previous result. The solid histograms use the same indices and similar priors as K03 or G05, while the 
dashed histograms use our set of indices and Monte-Carlo priors. 



Properties of LRGs in the SDSS 17 



G05 in to fairly close agreem ent. Had we adopted the lKroupal 
(|200ll ) or IChabrieJ l|2003l ) IMF, the peak in the measured 
stellar mass distribution would be smaller by about 0.2 dex, 
or 2 X 10" Mq. 

When using < 7 < 2, we also find that our best 
fitting models are more metal-rich compared with those of 
G05 by ~ 0.05 dex, with a scatter of ~ 0.15 dex. Our models, 
however, are comparatively metal-poor by ~ 0.05 dex when 
using < 7 < 1. Here we see the age-metallicity degeneracy 
in effect: younger ages correlate with higher metallicities, 
and vice versa. 

Sources of discrepancy in addition to the difference in 
priors include the rejection by K03 of models for which the 
r — i colour is redder than observed, the difference in the pop- 
ulation synthesis techniques adopted to generate the models 
and the consequent difference in the abilities of the models 
to reproduce the absorption line features in the spectra. We 
include models with attenuations down to —0.1, following 
the method of G05. We also note that we are using a differ- 
ent extinction curve to that adopted by K03 and G05, who 
both use the single power law (oc A~°'^) of Calzetti et al. 
l|l994 ). The derived Az values we find are sufficiently small 
(see Section [3. 4|l to have a negligible effect on the mass es- 
timates. 

One way to examine the effect of these different choices 
is to examine how well the two models reproduce t he ab - 
sorption indices. In figure 18 of iBruzual fc CharlotI (|2003l ). 
the H/3, USa, H7A, D„(4000), [MgFe]', [MgiFe] and [MgaFe] 
index strengths are used to fit high signal-to-noise galaxies 
from the SDSS Early Data Release (EDR). A comparison of 
this plot with our Figure [3] shows general consistency in the 
systematic offsets of many in dices, despite the broader range 
in galaxy types modelled bv lBruzual fc CharlotI (|2003l ) . For 
example, both models fit the Rf3, D„(4000) and [MgFe]' 
well, whilst underestimating the strengths of the CNi, CN2, 
Ca4455, C24668, Mg;, and NaD indices and overestimat- 
ing Fe4383, Fe5015, Fe5270 and Fe5335. This behaviour 
is encouraging as it allows us to check the consistency of 
both models. There are some indices that are not consis- 
tent, however. For example, we tend to overestimate the 
stre ngths of the two ma g nesiu m indices, Mgi and Mg2 , while 
the lBruzual fc CharlotI (|2003l ) models underestimate them. 
Other indices that show differences are Ca4227, G4300, 
Fe4531, Fe5406, TiOi, Ti02 and U5a- These differences are 
likely in part due to the differences in the galaxy populations 
modelled. 

The results of Section|4]4]highlight an important consid- 
eration when setting up this type of analysis, namely what 
range should 7 take in the prior? We have shown here that 
doubling the upper limit of 7 from 1 to 2 decreases the 
formation redshift considerably from 2 ~ 2.0 to 2 ~ 1.0, 
making the derived ages younger by around 2 Gyr. This 
anti-correlation between 7 and time since formation shown 
in Figure [5] is an example of the degeneracy between for- 
mation epoch and formation time-scale. So an old, slowly- 
forming galaxy would appear similar to a young galaxy 
that fo rmed over a relatively short time-scale. iThomas et al.l 
l|2005l ) claim that this degeneracy could be lifted by relating 
observationally derived a/Fe ratios with their time-scales. 
Their simulations predict that an early-type galaxy of mass 
M./M0 ~ 10"■^ should have formed on a time-scale of 
roughly 260 Myr. In their study, gaussians are used as a toy 



model of star formation histories, whereas we use decaying 
exponentials. Approximating their formation time-scale as 
the inverse of 7 gives 7 ~ 3.8 Gyr~^. This suggests that 
perhaps an even larger range of 7 is required to cover the 
full range of physical star formation histories. It appears, 
though, that the range we use works well for most of the 
galaxies. The distribution of 7 values in the upper panel of 
Figure [9] shows a gradual decline towards the upper limit 
for the < 7 < 2 fits. This is in contrast to a crowding 
of 7 values towards the upper limit for the < 7 < 1 fits. 
This suggests that, while some galaxies may be better fit by 
models with 7 > 2, the range < 7 < 2 well describes the 
bulk of the LRG population. 

We also check the consistency of our results by com- 
paring the trends shown by the 2D probability distributions 
of Figure [6] with previous studies. Figure 12 of G05 shows 
similar distributions of stellar metallicity and age for early- 
type galaxies in their sample. We are interested in the upper 
mass bin, where log(M/M0) > 11. They find that the high 
mass galaxies occupy a region of parameter space centred on 
an age of about 8 Gyr and metallicity of [Z/H] = 0.1 — 0.2. 
This is roughly consistent with the age-metallicity plot in 
Figure |6l which covers a similar region. Both plots also show 
a similar age metallicity-degeneracy trend, with a slope of 
AlogZ/Alogt = —1.05 dex yr"'^ from our plot, compared 
with the value of —0.75 dex yr"'^ found by G05. There are, 
however, some differences. We find ages smaller than theirs 
by about 1-2 Gyr and slightly higher metallicities. Our Fig- 
ure [12] demonstrates that these differences are due primarily 
to our different choice of indices to fit and Monte-Carlo pri- 
ors in the construction of the library of template spectra. 

Only a small percentage of the spectra are better fit 
by including a burst on top of the continuos star formation. 
The probability that a galaxy experienced a burst within 
the past 6 Gyr of its life is about 3 per cent, and only about 
0.25 per cent within the past 2 Gyr. This may in part be 
due to our selection criterion of avoiding spectra with emis- 
sion lines, but it is consistent with the small percentage of 
LRGs showing evidence of recent star formation found by 
[Roscboom fc et, al. ( 2006 ). For those galaxies better fit by 
the addition of a burst, the burst contributes 20-40 per cent 
of the total stellar mass of the galaxy. 

There have been a number of studies that show a dis- 
tinct trend between the properties of early -type galaxies and 
their velocity dispersions. Caldwell et al.l (|200S ') derive the 
ages and metallicities of 175 nearby, early-type galaxies us- 
ing a combination of index-index grids. They find a cor- 
relation between age and a, with a 230 kms"'^ (the up- 
per limit on their observations) galaxy expected to have 
an age between 4 and 14 Gyr and a metallicity range of 
—0.2 < [Z/H] < 0.3. Both of these results are broadly con- 
sistent with our findings, and suggest that we are studying 
too narrow a range in velocity dispersion or too homoge- 
neous a population to detect any significant trend. 

The evolution in the median star formation rate we in- 
fer for the LRGs agrees well with the semi-ana lytic galaxy 
formation predictions of iDe Lucia et al.1 (|2006D . We agree 
less well on the formation times. Comparing the times at 
which either 50 per cent or 80 per cent of the stars in the 
LRGs were formed gives times shorter by about 2 Gyr from 
the predictions of the semi-analytic models. We also find 
metallicities about 0.15-0.2 dex larger than predicted. 



18 Tom Barber, Avery Meiksin and Tara Murphy 



In an attempt to understand the origin of the metal- 
Hcity distribution of the galaxies, we introduce a stochastic 
star-formation quenching model for the metallicity evolu- 
tion. In this model, metallicity evolves either as a closed 
box or an accreting box, up until a time when star forma- 
tion abruptly ceases. The model is motivated by quenching 
scenarios invoked to resolve the 'anti-hierarchical' formation 
problem, such as the shutting-off of gas cooling by AGN ac- 
tivity. We find that while no model provides a perfect fit 
to the data, the underlying measured distribution is well- 
described by both the closed box and accreting box versions, 
provided a 'cuspy' burst of star-formation is adopted. The 
characteristic time-scale for quenching the star formation is 
around 10* yr, comparable to the estimated lifetimes of lu- 
minous QSOs. The required effective yield values, however, 
are quite different. For the closed box model, a yield of about 
0.26 is required. Such a high yield requires the enrichment 
to be dominated by very massive stars. Such a possibility, 
however, should perhaps not be discarded and is being con- 
sidered in a newer set of semi-analytic models (Lacey 2006, 
in preparation). For the accreting box, a more moderate 
yield of about 0.10 is required. This is still about twice the 
yield expected for the IChabrieij (|2003l ) IMF for early stars, 
and again suggests an IMF more strongly weighted towards 
higher masses. The dominant role Type II supernovae would 
play in determining the abundances for such IMFs is con- 
sistent with the Q-enhanced abundances inferred from the 
absorption indices of LRGs. 

The quenching model also predicts an anti-correlation 
between the ages of the stars today and their metallicities, as 
early quenching will result in older stars with smaller metal- 
licities for a given formation time. This provides a physical 
basis to the observed anti-correlation, which may thus not 
soley be an artefact of the age- metallicity degeneracy affect- 
ing population synthesis modelling. The relation, however, 
will also depend on the range of formation times and de- 
cay times for the star formation. Such effects are perhaps 
best accounted for in semi-analytic galaxy formation models. 
Here we suggest that the distribution of metallicity and the 
anti-correlation with age should provide further constraints 
on such models. 
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